Solution to Big-Bang Nucleosynthesis in Hybrid Axion Dark Matter Model 



O 
(N 

D 

in 



6 



> 
o 

o 



Motohiko Kusakabc^B A.B. Balantekin^'^ltl Toshitaka Kajino^'^jll and Y. Pchlivan^'H 
^Institute for Cosmic Ray Research, University of Tokyo, Kashiwa, Chiha 277-8582, Japan 
^Department of Physics, University of Wisconsin, Madison, WI 53706, USA 
'^National Astronomical Observatory of Japan 2-21-1 Osawa, Mitaka, Tokyo, 181-8588, Japan 
Department of Astronomy, Graduate School of Science, 
University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan and 
Mimar Sinan Fine Arts University, Besiktas, Istanbul 34349, Turkey 
(Dated: February 28, 2012) 

Following a recent suggestion of axion cooling of photons between the nucleosynthesis and recom- 
bination epochs in the Early Universe, we investigate a hybrid model with both axions and relic 
supersymmetric particles. In this model we demonstrate that the ^Li abundance can be consistent 
with observations without destroying the important concordance of deuterium abundance. 
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Introduction - Low-metallicity halo stars exhibit a 
plateau of ^Li abundance, indicating the primordial ori- 
gin of ^Li However, the amount of ^Li needed to be 
consistent with the cosmic microwave background obser- 
vations [3] is significantly more than ""Li observed in old 
halo stars (Even though ''Li can be both produced 
and destroyed in stars, old halo dwarf stars are expected 
to have gone through little nuclear processing). Recent 
improvements in the observational and experimental data 
seem to make the discrepancy worse 0, Mi ■ One possible 
solution is to invoke either nuclear physics hitherto ex- 
cluded from the Big-Bang Nucleosynthesis (BBN) d, Q 
or new physics such as variations of fundamental cou- 
plings Ip, |9l| , and particles not included in the Standard 
Model |lCll442l | . Effects of massive neutra l relic particles 
on BBN were extensively studied 
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More recently an alternative solution to the lithium 
puzzle was proposed. Since the last photon scattering oc- 
curs after the end of the nucleosynthesis, one can search 
for a mechanism for the cooling of photons before they 
decouple. It was suggested that dark matter axions could 
form a Bose-Einstein condensate (BEG) [43|, i4J] . Such 
a condensate would cool the photons between the end 
of BBN and epoch of photon decoupling, reducing the 
baryon-to-photon ratio WMAP infers, as compared to 
its BBN value [i^ . An alternative mechanism for such a 
cooling is resonant oscillations between photons and light 
abelian gauge bosons in the hidden sector ji^ . There are 
two prima facie problems with the axion BEC-photon 
cooling hypothesis: it overpredicts primordial deuterium 
(D) abundance as well as the effective number of neu- 
trinos. Even though D is easy to destroy, one does not 
expect the sum of abundances of D and '^He to change 
significantly in the course of cosmic evolution [47| . Hence 
it is important to find a parameter region in which pre- 
dicted abundances of D and ^Li are consistent with ob- 
servations. In this letter we demonstrate the existence of 
such a parameter region using a model with axions and 
massive relic particles. 



The Hybrid Model - We carried out BBN net- 
work calculations using Kawano's code (H^, [5l[ by in- 
cluding Sarkar's correction for ^He abundances [H^ . 
JINA REACLIB Database VI. ^ is used for light nu- 
clear {A < 10) reaction rates including uncertainties to- 
gether with data S^-l^l- Adopted neutron lifetime is 
878.5 ± O.Tstat ± O.Ssys s [13 based on improved mea- 
surements [4^. Taking into account the uncertainties in 
the rates of twelve important reactions in BBN |5l|, we 
employ regions of 95% C. L. in our calculations. 

We compare our results with the abundance con- 
straints from observations. For the primordial D abun- 
dance, the mean value estimated from Lyman-a absorp- 
tion systems in the foreground of high redshift quasi- 
stellar objects is log(D/H) = -4.55 ± 0.03 We 
adopt this value together with a 2a uncertainty, i.e., 
2.45 X 10-5 < D/H < 3.24 x 10"^ ^jje abun- 
dance measurements in Galactic HII regions through the 
8.665 GHz hyperfine transition of ■^He"'' yield a value of 
3Hc/H=(1.9 ± 0.6) X lO-'^ Although the constraint 
should be rather weak considering its uncertainty, we 
take a 2(t upper limit and adopt "^He/H < 3.1 x IQ-^. 

We also utilize a limit on the sum of primordial abun- 
dances of D and '^He taken from an abundance for the 
protosolar cloud determined from observations of solar 
wind, i.e., (D-h3He)/H=(3.6 ± 0.5) x 10"^ This 
abundance can be regarded as constant at least within 
the standard cosmology since it is not affected by stellar 
activities significantly despite an effect of D burning into 
^He via ^H(p, 7)^He would exist [47l |. 

For the primordial *He abundance, we adopt two dif- 
ferent constraints from recent reports: Yp = 0.2565 ± 
0.0051 [eH and Yp = 0.2561 ± 0.0108 [H both of which 
are derived from observations of metal-poor extragalactic 
HII regions. Adding 2a uncertainties leads to 0.2463 < 
Yp < 0.2667 [m and 0.2345 <Yp< 0.2777 [13 • 

®Li plateau of metal-poor halo stars (MPHSs), yields 
the upper limit of '^Li/H= (7.1 ±0.7) x lO^^^ Adding 
a 2a uncertainty, we adopt ^Li/H < 8.5 x 10"^^. 
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FIG. 1. Abundances of ''He (mass fraction), D, ^He, ^Li and 
®Li (number ratio relative to H) as a function of the baryon-to- 
photon ratio rj or the baryon energy density parameter ^Ish? 
of the universe. The thick dashed curves are for SBBN. The 
thin dashed curves around them show the regions of 95% C. 
L. in accordance with the nuclear reaction rate uncertainties. 
The boxes correspond to the adopted abundance constraints 
on the SBBN model. The vertical stripes represent the 2a 
limits on Qsh^ or t] for the SBBN model (taken from the 
constraint by WMAP ^ and labeled as WMAP7) and for 
the axion BEG model (labeled as axion). The solid curves 
are the results obtained with the long-lived decaying particle 
model with parameters fixed to {tx, (^x)~{W^ s, 2 x 10"'" 
GeV) (see text). 



For the ^Li abundance, we adopt the limits 
log(^Li/H)= -12 + (2.199 ± 0.086) (with 95% C. L.) de- 
rived from recent observations of MPHSs in the 3D non- 
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local thermal equilibrium model [63|, i.e. 1.06 x 10 
^Li/H < 2.35 X 10-1°. 

Figure [1] shows the abundances of ^He (^p; mass frac- 
tion), D, ^He, ^Li and ^Li (number ratio relative to H) as 
a function of the baryon-to-photon ratio 77 or the baryon 
energy density parameter Q,Bh^ of the universe, where h 
is the Hubble constant in units of 100 km/s/Mpc. The 



thick dashed curves are the results of the standard BBN 
(SBBN) with a neutron lifetime of 878.5 ± 0.8 s. Thin 
dashed curves around them show regions of 95% C. L. 
from uncertainties in the nuclear reaction rates. The 
boxes represent adopted abundance constraints as sum- 
marized above. The vertical stripes correspond to the 2a 
limits on ris/i^ or 77. The values provided by WMAP 
(labeled WMAP7) are 



(6.225l°ji4)xl0- 
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(1) 



Values predicted by the BEC model (labeled axion) are 
smaller by a factor of (2/3)3/^ at the BBN epoch 



(4.593l°:f,^)xl0-i°. (2) 



It can be seen that the adoption of the 77 value from 
WMAP leads to a ''Li abundance calculated in the ax- 
ion BEC model, which is in reasonable agreement with 
the observations. However, we lose the important consis- 
tency in D abundance. Ref. [45| noted that astronom- 



ical measurements of primordial D abundance can have 
a significant uncertainty as well as a possibility that D 
is burned by nonstandard stellar processes. Even if their 
assumption were true, stellar processes are not expected 
to change the sum of D and "^He abundances 47[. As 
seen in Fig. [U the constraint on (D-|-'^He)/H abundance 
seems to exclude the original axion BEC model. Ulti- 
mately, this model is viable only when the abundance of 
(D+'^He)/H is reduced through some exotic processes. 

It is known that nonthermal photons can be generated 
through electromagnetic energy injections by the radia- 
tive decay of long-lived particles after the BBN epoch [ll|, 
[Tij . These nonthermal photons can photodisintegrate 
background light elements ll( 
adopt the method of Ref. [251 




We 



to calculate the nonther- 
mal nucleosynthesis, where we incorporated new thermal 
reaction rates as described above. In addition, we adopt 
updated reaction rates of "He photodisintegration [26[ 
derived from the cross section data using precise mea- 
surements with laser-Compton photons [6J,|63. Effects 
of electromagnetic energy injection depend on two pa- 
rameters. One is C,x — k^xl"^-^^!^ where (n^/n°) is 
the number ratio of the decaying particle X and the back- 
ground radiation before the decay of X, and -E-yo is the 
energy of photon emitted at the radiative decay. The 
other is tx , the lifetime of the X particle. 

Figure [2] shows the parameter space in our hybrid 
model. For "He, we adogt the conservative constraint 
with larger uncertainty j62l |. We also show the contour 
for ^Li production at the observed level, i.e., ^Li/H= 
7.1 x 10^1^. This figure shows the result of nonthermal 
nucleosynthesis induced by the radiative decay of long- 
lived particles with the 77 value of the axion BEC model 
[Eq. Except for D and ''Li, contours are similar to 

those represented in Ref |25| where BBN epoch 77 value 
is assumed to be the same as WMAP i] value [Eq. ([T])]. 
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FIG. 2. Parameter space of the hybrid model (rx, Cx) for the 
value of 77 = 4.6 x 10"^" provided by the axion BEG model. 
The contours identify the regions where the nuclei are over- 
produced or underproduced ("over" and "low", respectively) 
with respect to the adopted abundance constraints. "^He mass 
fraction (red line) and ^He/H (black lines), D/H (green solid 
and dashed lines for upper and lower limits, respectively), 
and ''Li/H (blue line) number ratios are shown. The orange 
line is the contour of ^Li/H=7.1 x 10~^^. In the gray-colored 
region all abundances are within the limits of observational 
constraints. 



In the very small colored region, calculated primordial 
abundances of all nuclides including D and ^Li are simul- 
taneously in ranges of adopted observational constraints. 
We conclude that the present model eliminates the main 
drawback of the original axion BEC model by reducing 
primordial D abundance via ^ll{'j,n)^ll reaction, where 
7's are nonthermal photons. We note that the decaying 
particle model with the WMAP 77 value cannot resolve 
the ^Li problem by itself (HIH]. 

The effect of the radiative decay on other elemental 
abundances is not significant except for ''Li. Since ener- 
getic photons produced quickly collide with background 
photons and create e~^e~ pairs, nonthermal photon spec- 
tra developed by the decay has a cutoff energy Eq = 
mg/22r where me is the electron mass The decay 
at earlier (hotter) universe then triggers nonthermal pho- 
tons with lower cutoff energies. The threshold energies 
of ^Be and D photodisintegration, ^Be-l-7 — >^He+'*Hc 
and B+j n + p are 1.5866 MeV and 2.2246 MeV, 
respectively. These two nuclei are very fragile against 
photodisintegration. When the decay occurs early at rel- 
atively high Tg = r/(10^ K) > 10-2 which corresponds 
to Tx ^ 10^ s, nonthermal photon spectra contain pho- 
tons to dissociate ^Be and D, while keeping other nuclides 
intact. The gray region indicates parameters which re- 
sult in moderate destruction of D which is overproduced 
in the original BEC model because of low rj. Above that 
region, D is destroyed too much by photodisintegration. 




FIG. 3. Mass fractions of H and *He {Xp and Yp, respec- 
tively) and number ratios of other nuclides relative to H as 
a function of Tg. Solid lines show the abundances calculated 
in the hybrid model with the parameters (rx, Cx) = (10^ s, 
2 X 10"^" GeV) which correspond to the point indicated with 
a star in Fig. [S] The dashed lines show the SBBN prediction. 



while below it D abundance is too high. 

We next present the results of a BBN calculation in our 
hybrid model with a fixed set of parameters given by (rx , 
(^) = (106 s, 2 X 10^1" GeV) and the rj value given in Eq. 
([2| . This corresponds to the point indicated with a star in 
Fig. ^ and yields the required D and ^Li abundances as 
seen in Figure [31 This figure shows H and ''^He mass frac- 
tions (denoted by Xp and Yp), and n, D, T, '^He, ^Li, '^Li 
and ''Be number ratios relative to H as a function of the 
temperature. The abundances calculated using the axion 
BEC -I- long-lived decaying particle model with the pa- 
rameters (t_x, (x)={W^ s, 2x 10"^° GeV) and the 77 value 
provided by BEC model [Eq. ([2])] are shown in solid lines 
whereas the SBBN prediction is plotted by dashed lines. 
The small difference at Tg > 0.06 observed between solid 
and dashed lines is caused by difference between initial 
77 values. At 0.06 ^ Tg > 7 x lO""^ (corresponding to the 
cosmic time of t 5 x 10^^-4 x 10^ s), effects of ^11(7, ?i)^H 
are seen in the decrease of D and the increase of n abun- 
dances. We find a slight decrease in ''Be abundance. 
This is caused through reactions ^Be(7,'^He)^He (thresh- 
old energy: 1.5866 MeV), ^Be(7,p)6Li (5.6858 MeV), and 
'^Be(7, 2p??,)4He (9.3047 MeV). The second reaction in- 
creases the ^Li abundance. Finally, at Tg ^ 7 x 10~^, 
where the abundance of long-lived X particle is already 
less than 3 % of the initial abundance, effect of ''He pho- 
todisintegration is to increase "^H and n abimdances. 

In Fig. [1] solid lines show the results for the SBBN 
+ long-lived decaying particle model with the same pa- 
rameter values, i.e., {tx, Cx)=(10^ s, 2 x 10"^° GeV). 
Obviously the abundances of D and ^Li (produced partly 
as ^Be) are reduced, while that of ^Li is increased from 
those of SBBN. 
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Conclusions - We used a hybrid axion and massive 
relic particle model in which axions cool the photons and 
relic particles produce non-thermal photons to eliminate 
the high D abundance in the original axion BEC model. 
Our hybrid model also produces ^Li keeping ^Li abun- 
dance at the level of Population II Spite plateau. Our 
work thus demonstrates that the ^Li abundance can be 
consistent with observations without destroying the im- 
portant concordance of deuterium abundance. 
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